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Abstract 

We show that the evaluation of scattering ampHtudes can be formulated as a problem of multivariate polyno- 
mial division, with the components of the integration-momenta as indeterminates. We present a recurrence 
relation which, independently of the number of loops, leads to the multi-particle pole decomposition of 
the integrands of the scattering amplitudes. The recursive algorithm is based on the Weak Nullstellensatz 
Theorem and on the division modulo the Grobner basis associated to all possible multi-particle cuts. We 
apply it to dimensionally regulated one-loop amplitudes, recovering the well-known integrand-decomposition 
formula. Finally, we focus on the maximum-cut, defined as a system of on-shell conditions constraining the 
components of all the integration-momenta. By means of the Finiteness Theorem and of the Shape Lemma, 
we prove that the residue at the maximum-cut is parametrized by a number of coefficients equal to the 
number of solutions of the cut itself. 
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1. Introduction 

Scattering amplitudes in quantum field theories 
are analytic functions of the momenta of the par- 
ticles involved in the scattering process, and can 
be determined by their singularity structure. The 
multi-particle factorization properties of the ampli- 
tudes are exposed when propagating particles go on 
their mass-shell [H-Q- 
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The investigation of the residues at the singular 
points has been fundamental for discovering new 
relations fulfilled by scattering amplitudes. The 
BCFW recurrence relation [^], its link to the lead- 
ing singularity of one-loop amplitudes [3], and the 
OPP integrand-decomposition formula for one-loop 
integrals (5|] have shown the underlying simplicity 
beneath the rich mathematical structure of quan- 
tum field theory. Moreover they have become effi- 
cient techniques leading to quantitative predictions 
at the next-to-leading order in perturbation theory 

The integrand reduction methods [B*] allow to de- 



compose one-loop amplitudes in terms of Master 
Integrals (MPs) without performing the loop inte- 
gration, and are based on the multi-particle pole ex- 
pansion of the integrand. The expansion is equiva- 
lent to the decomposition of the numerator in terms 
of (a combination of) products of denominators, 
with polynomial coefficients. In the context of an 
integrand-reduction, any integration is replaced by 
polynomial fitting. 

The first extension of the integrand reduction 
method beyond one- loop was proposed in [14], and 
it was used to reproduce the results of two-loop 5- 
point planar and non-planar amplitudes in Af — A 
SYM |ll,[3. A key point of the higher-loop exten- 
sion is the proper parametrization of the residues 
at the multi-particle poles. Each residue is a mul- 
tivariate polynomial in the irreducible scalar prod- 
ucts (ISP's) among the loop momenta and either ex- 
ternal momenta or polarization vectors constructed 
out of them. ISP's cannot be expressed in terms of 
denominators, thus any monomial formed by ISP's 
is the numerator of a potential MI which may ap- 
pear in the final result. Hence, a systematic classi- 
fication of the polynomial structures of the residues 
is mandatory. In [l3 |. the residues have been ob- 
tained by relating the ISPs to monomials in the 
components of the loop momenta expressed in a ba- 
sis chosen according to the topology of the on-shell 
diagram. 

Badger, Frellesvig and Zhang [17,] combined on- 
shell conditions with Gram-identities [lH to limit 
the number of monomials appearing in the residues. 
This technique was applied to the integrand de- 
composition of two-loop 4-point planar and non- 
planar diagrams in supersymmetric as well as non- 
supersymmetric YM theories. 

In this work, we show that the shape of the 
residues is uniquely determined by the on-shell con- 
ditions alone, without any additional constraint. 
We derive a simple integrand recurrence relation 
that generates the required multi-particle pole 
decomposition for arbitrary amplitudes, indepen- 
dently of the number of loops. 

The algorithm treats the numerator and the de- 
nominators of any Feynman integrand, as multi- 
variate polynomials in the components of the loop 
variables. The properties of multivariate polyno- 
mials have been extensively studied in the math- 
ematical literature, see e.g. (l9l - (25} . The method 
uses both the weak Nullstellensatz theorem and the 
multivariate polynomial division modulo appropri- 
ate Grobner basis [T^. In the context of the in- 



tegrand reduction, univariate polynomial division 
has been already introduced in [26| to improve the 
decomposition of one-loop scattering amplitudes. 

The algorithm, which is described in Section II, 
relies on general properties of the loop integrand: 

• When the number n of denominators is larger 
than the total number of the components of 
the loop momenta, the weak Nullstellensatz 
theorem yields the trivial reduction of an n- 
denominator integrand in terms of integrands 
with (n — 1) denominators. 

• When n is equal or less than the total number 
of components of the loop momenta, we divide 
the numerator modulo the Grobner basis of the 
n-ple cut, namely modulo a set of polynomi- 
als vanishing on the same on-shell solutions as 
the cut denominators. The remainder of the 
division is the residue of the n-ple cut. The 
quotients generate integrands with (n — 1) de- 
nominators which should undergo the same de- 
composition. 

• By iterating this procedure, we extract the 
polynomial forms of all residues. The algo- 
rithm will stop when all cuts are exhausted, 
and no denominator is left, leaving us with the 
integrand reduction formula. 

In Section III we apply the algorithm to a generic 
one-loop integrand, reproducing the d-dimensional 
integrand decomposition formula 0, [27l - [29| . 

In Section IV we conclude by proving a theorem 
on the maximum- cuts, i.e. the cuts defined by the 
maximum number of on-shell conditions which can 
be simultaneously satisfied by the loop momenta. 
The on-shell conditions of a maximum cut lead 
to a zero-dimensional system. The Finiteness 
Theorem and the Shape Lemma ensure that the 
residue at the maximum-cut is parametrized by Us 
coefficients, where is the number of solutions 
of the multiple cut-conditions. This guarantees 
that the corresponding residue can always be 
reconstructed by evaluating the numerator at the 
solutions of the cut. 

During the completion of this work, Zhang has 
presented an algorithm [sO] embedding the ideas 
presented in fl7t within more general techniques of 
algebraic geometry, among which the division mod- 
ulo Grobner basis is used as well. 
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2. Multivariate polynomial division 

In what follows, we assume 4-diniensional loop- 
momenta. Extensions to higher-dimensional cases, 
according to the chosen dimensional regularization 
scheme, can be treated analogously - as we will show 
when discussing the one-loop integrand reduction. 

The integrand reduction methods 0, [H [H, [26l - 
H, 31 "34 1 recast the problem of computing £-\oop 
amplitudes with n denominators as the reconstruc- 
tion of integrand functions of the type 



Ai(gi, ■■■,qe)- ■■D^^{qi, ...,qi) 



, (1) 



where qi, . . . ,qi are integration momenta. The 
generic propagator can be written as follows: 



A = a, + p,j - m\ , a, £ {0, ±1} . (2) 

The numerator Afi^...i^ and any of the denominators 
Di are polynomial in the components of the loop 
momenta, say z = (zi, . . . , Z4£), i.e. 



J. 



Mi-»,.(z) 
Ai(z)--- A„(z) 



(3) 



Let us consider the ideal generated by the n de- 
nominators in Eq. (jSj , 

\Ji\ ■■■in {-^il 5 ' ' " 5 ) 

= |iZ^-(^)A„(z):/i«(z)eP[z]|, 

where P[z] is the set of polynomials in z. The com- 
mon zeros of the elements of i7ii---i„ are exactly the 
common zeros of the denominators. 

The multi-pole decomposition of Eq. ^ is ex- 
plicitly achieved by performing multivariate poly- 
nomial division, yielding an expression of Afi-^...i^ in 
terms of denominators and residues. 
We construct a Grobner basis (see Ch. 2 of 
(ioj'). generating the ideal J7ii...i„ with respect to a 
chosen monomial order. 



^^■■■t„ = • • ■ ,gm(z)} 



(4) 



Unless otherwise indicated, we will assume lexico- 
graphic order. 

In this formalism, the n-ple cut-conditions = 
. . . = Di„ = 0, are equivalent to 171 = . . . = = 0. 
The number m of elements of the Grobner basis is 



the cardinality of the basis. In general, m is dif- 
ferent from n. We then consider the multivariate 
division of 7Vii...i„ modulo ^ii - i„ (see Ch. 2, Thm. 
3 of M), 



(5) 



where Ti-^...i^ = Qi(z).9i(z) is a compact nota- 

tion for the sum of the products of the quotients Qi 
and the divisors gi. The polynomial A^^...^^ is the 
remainder of the division. Since Gi-^...i^ is a Grobner 
basis, the remainder is uniquely determined once 
the monomial order is fixed. 

The term Ti-^...i^ belongs to the ideal Ji^...i^, thus 
it can be expressed in terms of denominators, as 

n 

K=l 

The explicit form of A/ij...i^_j,;^^j...i„ can be found 
by expressing the elements of the Grobner basis in 
terms of the denominators. 

2.1. Reducibility criterion. 

An integrand Xi-^...i^ is said to be reducible if it 
can be written in terms of lower-point integrands: 
that happens when the numerator can be written in 
terms of denominators. The concept of reducibility 
can be formalized in algebraic geometry. Indeed a 
direct consequence of Eqs. ([5|) and (O is the follow- 
ing 

Proposition 2.1. The integrand is re- 

ducible iff the remainder of the division modulo a 
Grobner basis vanishes, i.e. iff Mi-^...i^ G Ji^-.-i^- 

Proposition 12.11 allows to prove 

Proposition 2.2. An integrand 1,1, is re- 
ducible if the cut (ii,...,i„) leads to a system of 
equations with no solution. 

Proof. In this case, the system is over-constrained. 
The n propagators cannot vanish simultaneously, 
i.e. 

A,(z) = •■• = A„(z) =0 (7) 

has no solution. Therefore, according to the weak 
Nullstellensatz theorem (Thm. 1, Ch. 4 of 20]), 



K=l 



(8) 
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for some w„ G -P[z]. Irrespective of the monomial 
order, a (reduced) Grobner basis is Q = {gi} — {!}. 
Eq. (O becomes 



(z) =M,...,„(z) X 1 e Ji^...,,^ , (9) 



thus is reducible. 



□ 



2.2. Integrand Recursion Formula 

After substituting Eqs. ^ and ([6]) in Eq. we 
get a non-homogeneous recurrence relation for the 
n-denominator integrand, 



■A, 



(10) 



According to Eq. (fTO|) . li^-.-i^ is expressed in terms 
of integrands, Iii...i„_ii„^.ii„, with (n — 1) denomi- 
nators, are obtained from 2^^...^,^ by 
pinching the i^-th denominator. The numerator of 
the non- homogeneous term is the remainder A^j^...^^ 
of the division ([5]) . By construction, it contains only 
irreducible monomials with respect to Gij^---i„, thus 
it is identified with the residue at the cut (zi . . . in)- 

The integrands Iii...i„_ii„+i...i„ can be de- 
composed repeating the procedure described in 
Eqs. ([n])-®. In this case the polynomial division 
of has to be performed modulo the 

Grobner basis of the ideal J7ii...i^_ii^^i...i„ , gener- 
ated by the corresponding (n — 1) denominators. 

The complete multi-pole decomposition of the in- 
tegrand 2,1 is achieved by successive iterations 
of Eqs. ©-(IS]). Like an Erathostene's sieve, the re- 
cursive application of Eqs. ([5]) and (ITO)) extracts the 
unique structures of the remainders A's. The pro- 
cedure naturally stops when all cuts are exhausted, 
and no denominator is left, leaving us with the in- 
tegrand reduction formula. 

If all quotients of the last divisions vanish, the 
integrand is cut-constructible, i.e. it can be deter- 
mined by sampling the numerator on the solutions 
of the cuts. If the quotients do not vanish, they 
give rise to non-cut-constructihle terms, i.e. terms 
vanishing at every multi-pole. They can be recon- 
structed by sampling the numerator away from the 
cuts. Non-cut-constructible terms may occur in 
non-renormalizable theories, where the rank of the 
numerator is higher than the number of denomina- 
tors [Hi. 



The Proposition 12.21 and the recurrence rela- 
tion (jlOl) are the two mathematical properties un- 
derlying the integrand decomposition of any scat- 
tering amplitudes. The polynomial form of each 



residue is univocally derived from the division mod- 
ulo the Grobner basis of the corresponding cut. 

3. One-loop integrand decomposition 

In this section we decompose an n-point inte- 
grand Xq...(„_2) of rank-n with n > 5, using the 
procedure described in Section [51 The reduction of 
higher-rank and/or lower-point integrands proceeds 
along the same lines. 

In d-dimensions, the generic n-point one-loop in- 
tegrand reads as follows: 



O---(n-l) 



M)...(„_i)(g,At^) 
-Do(<7, M^) • • •^n-i(<7, M^) 



(11) 



We closely follow the notation of [26|, |35| . Objects 
living in d = 4 — 2e are denoted by a bar, e.g. 
I = ^ + /i and = — 

For later convenience, for each Xii-.-i^. we define 
a basis f — {ei}i=i^,,,^i. 

If fc > 5, then Ci = ki, where ki are four external 
momenta. 

If fc < 5, then Ci are chosen to fulfill the following 
relations: 



e? = = 



ei • 62 = 1 , 



63 • 64 



-(I-44). (12) 



In terms of ■ the loop momentum can be 
decomposed as, 

= ~Pi^ + xi Ci + X2 6^ -I- X3 ei^ -|- X4 6^ . (13) 

Accordingly, each numerator J\fi-^...i^, can be treated 
as a rank-fc polynomial in z = {xi,X2, X3,X4, ^'^), 



J\f^ ^ — \ .,J4 
jGJ(fc) 



^ aj z{- z^' zf zi- z^" , (14) 



with J(fc) = {/ = (ji, . . . , js) : ji -I- i2 + ja + 3i + 
2j5 < fc}. 

Step 1. When n > 5, the Proposition 12.21 guar- 
antees that A/o...n~i is reducible, and, by iteration, 
it can be written as a linear combination of 5-point 
integrands 

Step 2. The numerator of each is a rank-5 

polynomial in z, cfr. Eq. (fH)) . We define the ideal 
J^i-^...ij^, and compute the Grobner basis Qi-^...i^ — 



4 



. . . , (75), which is found to have a remarkably 
simple form: 

g.i{z) ^ Ci + Zi , (i = 1, . . . , 5) . (15) 

We observe that each gi depends linearly on the i-th 
component of z. 

The division of A/ij-.-ij; modulo Gii---ii,, see Eq.(IS]), 
gives a constant remainder, 



Ai^...ig — Co 

The term Ti-^...i^ in Eq. ^ is, 



(16) 



where A/ii...i„_ii^^i...i5 are the numerators of the 4- 
point integrands, Xi-^...i^_-^i^_^-^...i^, obtained by re- 
moving the i^-th denominator. 

Step 3. For each 1^^...;^ , the numerator A/ii...i4 is 
a rank-4 polynomial in z. The Grobner basis Qi^...i^ 
of the ideal Ji^...i^ contains four elements. Divid- 
ing A/ii...i4 modulo Gii---ii, we obtain the remainder. 
The latter depends on /x^ and on the fourth compo- 
nent of the loop momentum q in the basis £:(*i---i4)^ 



2/^ I „ ^ 1 ,,2, 



+ Ai (C2 + €3X4 + fjL C4) . 



(17) 



The term F,- 
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A/ij...i^_ji^^4...i4 (z)_Di^ (z) , 



contains the numerators of 3-point integrands 
7. . 

S'fep The Grobner basis Gixi^i^ is formed 
by three elements, and is used to divide Ni^i^i^. 
The remainder Ai^i^i^ is polynomial in fi^ and in 
the third and fourth components of q in the basis 
£■(111213) 

^ililh = Co + CiXz + C2xl + C^xl 
+ C4X4 + C^x\ + Cqx\ 

+ ^^(07 + 08X3 + 09X4). (18) 

The term F^ji^ig generates the rank-2 numerators 
of the 2-point integrands T-i^i^^ '^iii^i and 1^2*3 ■ 

Step 5. The remainder of the division of A/i^ij by 
the two elements of Qi^i^ is: 

^^7,1 Zi 



= Co + cia:;2 + 02X3 + C3X4 
+ Cix\ + 05X3 + cqx\ + 07X2X3 

+ CgX2X4 + Cg^^ . 



(19) 



It is polynomial in /i^ and in the last three compo- 
nents of q in the basis ^('i*^) , xhe reducible term of 
the division, Fij^^, generates the rank-1 integrands, 
li^ , and . 

Step 6. The numerator of the 1-point integrands 
is linear in the components of the loop momentum 
in the basis 



4 

The only element of the Grobner basis Qi^ \s Di^, 
which is quadratic in z. Therefore the division mod- 
ulo Qi^, leads to a vanishing quotient, hence 



Ml = A,, 



(20) 



Step 7. Collecting all the remainders computed 
in the previous steps, we obtain the complete de- 
composition of To. ..,1-1 in terms of its multi-pole 
structure 



0---n-l — 




A,; 



D, 



■D, 



(21) 



Eq. ((2T|) reproduces the well-known one-loop d- 
dimensional integrand decomposition formula 0, 

We remark that the basis f ('i" *'--), defined in 
Eg.dT^ and used for decomposing the integration 
momentum q, depends only on the external mo- 
menta of diagram associate to the cut, eventually 
complemented by orthogonal elements. Therefore, 
£(ii---'ik) (.g^jj j-jg uged as well to decompose the inte- 
gration momenta of multi-loop diagrams (l^ . 



4. The Maximum- cut Theorem 

At £ loops, in four dimensions, we define a 
maximum- cut as a (4£)-ple cut 



D,. 



D, 







(22) 



which constrains completely the components of the 
loop momenta. In four dimensions this implies the 
presence of four constraints for each loop momenta. 
We assume that, in non-exceptional phase-space 
points, a maximum-cut has a finite number of 
solutions, each with multiplicity one. Under this 
assumption we have the following 
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Theorem 4.1 (Maximum cut). The residue at the 
maximum- cut is a polynomial paramatrized by Ug 
coefficients, which admits a univariate representa- 
tion of degree (n^ — 1)- 

Proof. Let us parametrize the propagators using M 
variables z = (2:1, ... , zu). In this parametrization, 
the solutions of the maximum-cut read, 



(4^ 



withi — 1, 



(23) 



Let Jii-.-i^i be the ideal generated by the on-shell 
denominators, Ji^...i^, = (Ai , ■ • ■ , A^) ■ 
According to the assumptions, the number rig of the 
solutions of (211) is finite, and each of them has mul- 
tiplicity one, therefore Ji^-.-is^^ is zero-dimensional 
[11 and radical, see Cor. 2.6, Ch. 4 of [21^. In 
this case, the Finiteness Theorem (Prop. 8, Ch. 
5 of [2^) ensures that the remainder of the di- 
vision of any polynomial modulo Ji-^^...i^^ can be 
parametrized exactly by coefficients. 

Moreover, up to a linear coordinate change, we 
can assume that all the solutions of the system 
have distinct first coordinate zi, i.e. z^^ ^ z'"'^ 

i ^ j. We observe that Ji^...i^^ and z\ are in the 
Shape Lemma position (Prop. 2.3 of i24j ) there- 
fore a Grobner basis for the lexicographic order 
zi < Z2 < ■ ■ ■ < Zn IS Qi^...i^, = {51,..., in 
the form 



51 (z) 
52(2) 

ff«(z) 



/i(^i) 

Z2 " f2{z\) 

Zit ~ fuizi) 



(24) 



The functions fi are univariate polynomials in zi. 
In particular /i is a rank-ris square-free polynomial 



It's 



(25) 



i.e. it does not exhibits repeated roots. The multi- 
variate division of A/ii...i4^ modulo Gii-.-iu leaves a 
remainder A^j^...^^^ which is a univariate polynomial 



in zi of degree (n^ — 1) [2£ 
Finiteness Theorem. 



in accordance with the 
□ 



The maximum-cut theorem ensures that the 
maximum-cut residue, at any loop, is completely 



^ Given an ideal JT, the radical of J7 is = {/ 6 P[z] 
3 s 6 N, e J}. J is radical iS J = y/j. 
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Figure 1: The on-shell diagrams in the picture are exam- 
ples of maximum-cuts. The first diagram in the left column 
represents the 5-ple cut of the 5-point one-loop dimension- 
ally regulated amplitude. All the other on-shell diagrams are 
considered in four dimensions. For each of them, the general 
structure of the residue A (according to the Shape Lemma) 
and the corresponding value of Ug are provided. 



determined by the distinct solutions of the cut- 
conditions. In particular it can be reconstructed 
by sampling the integrand on the solutions of the 
maximum cut itself. 

At one loop and in (4 — 2e)-dimensions, the 5- 
ple cuts are maximum-cuts. The remarkably sim- 
ple structure of the Grobner basis in Eq. (1X5]) is 
dictated by the maximum-cut theorem. Moreover 
in this case Ug = I, thus the residue in Eq. p6p is 
a constant. 

The structures of the residues of the maximum 
cut, together with the corresponding values of Us, 
for a set of one-, two-, and three- loop diagrams are 
collected in Figure [TJ 

The calculations of Sections [3] and p have been 
carried out using the package S@M [37| and the 
functions GroebnerBasis and PolynomialReduce 
of Mathematica, respectively needed for the gen- 
eration of the Grobner basis and the polynomial 
division. 

5. Conclusions 

We presented a new algebraic approach, where 
the evaluation of scattering amplitudes is addressed 
by using multivariate polynomial division, with 
the components of the loop-momenta as indetermi- 
nates. We found a recurrence relation to construct 
the integrand decomposition of arbitrary scattering 
amplitudes, independently of the number of loops. 
The recursive algorithm is based on the Weak NuU- 
stellensatz Theorem and on the division modulo 
the Grobner basis associated to all possible multi- 
particle cuts. Using this technique, we rederived 
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the well-known one-loop integrand decomposition 
formula. Finally, by means of the Finiteness Theo- 
rem and of the Shape Lemma, we proved that the 
residue at the maximum-cuts is parametrised ex- 
actly by a number of coefficients equal to the num- 
ber of solutions of the cut itself. 
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